Measurement of the resonant dfit molecular formation rate in solid HD 
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Measurements of muon catalyzed dt fusion (dfit —> 4 He + 
n + fi~) in solid HD have been performed. The theory 
describing the energy dependent resonant molecular forma- 
tion rate for the reaction /it + HD — > [(dfit)pee\* is com- 
pared to experimental results in a pure solid HD target. 
Constraints on the rates are inferred through the use of a 
Monte Carlo model developed specifically for the experiment. 
From the time-of-flight analysis of fusion events in 16 and 37 
/ig-cm -2 targets, an average formation rate consistent with 
0.897±(0.046) S iat ± (0.166) S yst times the theoretical predic- 
tion was obtained. 
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The idea of introducing a negative muon into a mixture 
of hydrogen isotopes to catalyze fusion has been known 
for over 50 years [Q. The muon will be quickly cap- 
tured by an atom, replacing the orbital electron. This 
muonic atom can then interact with an isotopic hydro- 
gen molecule to form a muonic molecular ion. The neg- 
ative muon is about 207 times the mass of an electron, 
thus the two nuclei in the molecule are bound together 
very tightly. The negative muon shields the repulsive 
Coulomb force between the two nuclei and because of 
their close proximity, fusion occurs and the muon can 
take part in another fusion cycle. The rate-limiting step 
in this process is muonic molecular (or molecular ion) for- 
mation. However, the probability of the muon sticking to 
the alpha particle after fusion occurs limits the average 
number of fusions to approximately 200 per muon. 

The theory describing the energy dependent resonant 
molecular formation rate for the reaction fit + HD — > 
[(dfit)pee\* has been developed in the last decade. The 
high rate is expected (|] to enhance muon catalyzed fu- 



sion in triple mixtures of H/D/T. In this letter we re- 
port the results of a muon catalyzed dt fusion experi- 
ment which confirm this rate. We have performed the 
first measurement using fit atoms incident on pure solid 
HD, measuring energy via a time of flight technique. By 
using a Monte Carlo simulation (which employs the the- 
oretical energy-dependent rates) we have compared the 
experimental results to the theory, thus determining the 
resonant molecular formation rate consistent with the fu- 
sion time spectra for two different HD target thicknesses. 

Figure [l] shows the theoretically predicted resonant 
molecular formation rate as a function of fit lab energy. 
This calculation, normalized to liquid hydrogen atomic 
density (LHD), was performed for isolated HD molecules 
at a temperature of 3K as had been done at higher tem- 
peratures H . The rate for HD is then predicted by multi- 
plying by the fractional deuteron concentration, Cd=0.5. 
The excess collisional energy of the system along with 
the released binding energy of the dfit molecular ion is 
absorbed by excitation of the entire six-body complex 
[(dfit)pee\. Because the dfit ion exists in a weakly bound 
state, resonance formation is possible without dissocia- 
tion. 

The experiment described here has the advantage of 
using spatially separated layers of different isotopic com- 
position in which interactions of the muonic system can 
be separated, providing the opportunity to measure the 
rate and energy dependence of resonant molecular for- 
mation using time of flight. It is similar to our previ- 
ous reported measurement Q using D 2 molecules. Other 
experimental efforts in liquid/gaseous triple mixtures at- 
tempt to measure this rate, however interpretation is dif- 
ficult because of the complex kinetics involved [5,6]. In 
the present experiment, solid layers were made by freez- 
ing hydrogen isotopic gases and mixtures onto two 51 fim 
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thick gold target foils maintained at 3K. The gas mix- 
tures were deposited onto the foils through a gas diffuser 
[7,8], where they could remain frozen for up to one week. 
For the HD target layer, commercial 96% HD gas was 
passed through a molecular sieve placed in a liquid ni- 
trogen bath to remove nonhydrogen impurities, and was 
subsequently monitored during target deposition with a 
quadrupole mass spectrometer to confirm noncquilibra- 
tion of the gas. It has been shown that at cryogenic tem- 
peratures the molecules do not begin to equilibrate for 
roughly 200 hours [pj . A nonhydrogen contamination fol- 
lowing target deposition was determined to be less than 
0.1% through use of a germanium detector which mea- 
sured x-rays from muon transfer to high-Z nuclei as well 
as gamma rays from nuclear capture. 

Figure || depicts the target and the muonic interac- 
tions which occur. Negative muons having a momentum 
near 27 MeV/c travel through one of the gold foils (a) 
into the 3.5 mg-cm~ 2 production layer (b) which consists 
of protium ^Eb) containing 0.1% tritium (T2). Here 
the muon is atomically captured by a proton within a 
time of 10~ n s (at LHD) and subsequently transfers to 
a triton due to the increased binding energy of fit. The 
neutral fit atoms, created with tens of eV energy, es- 
cape the production layer before thermalization due to 
the Ramsauer-Townsend mechanism which reduces their 
scattering cross section on H 2 . The moderation layer (c), 
consisting of 70 /ig-cm -2 D2, reduces the energy of the 
lit atoms to optimize resonant molecular formation on 
the downstream HD target (d) of thickness 16 /ig-cm~ 2 . 
The fit atoms traverse 17.9 ± 0.1 mm in vacuum to reach 
(d). For those at or near the resonant energy, molecular 
formation and fusion follow with high probability, pro- 
ducing a 3.5 MeV alpha particle and a 14 MeV neutron. 
The fit atoms which do not escape the moderation layer 
(c) can interact there with deuterium, and subsequent dt 
fusion can occur 

Alpha particles from fusion in (c) and (d) are separated 
by time of flight. The muonic atom emission time from 
the production layer, dominated by muon transfer from 
proton to triton, is of order 0.1 /is. Resonant formation 
rates are of order 10 9 s _1 , while dt fusion occurs in the 
molecular ion at a rate of 10 12 s — 1 . Since the two gold 
foils are separated by 17.9 mm, a fit atom which travels 
through the moderator layer reaches the downstream HD 
layer in a minimum time of 0.2 /is, and more typically 2 
[is at 1 eV. Thus, the time interval between detection of a 
muon entering the target and an alpha from fusion in the 
downstream HD reaction layer is dominated by the flight 
time of the fit atom. Two planar silicon detectors located 
perpendicular to the target foils, each with active area 
2000 mm 2 and nominal depletion depth 300 fim, detect 
alpha particles. A time cut of 1.5 to 6.0 /is, selects fu- 
sion events occurring in the HD layer only. The summed 
energy spectrum from both silicon detectors, normalized 
to the number of good incident muons (see for 



details), is shown in Fig. ||. Background was removed by 
subtracting normalized data for which no downstream 
reaction layer was present. 

A Monte Carlo (MC) program was used to compare 
theoretical resonant molecular formation rates to data. 
This was necessary as the time-of-flight spectra cannot 
be uniquely inverted due to geometrical effects and the 
energy loss of fit atoms in the HD reaction layer prior 
to molecular formation. The MC code |l(| simulates 
muonic processes occurring in the experiment (e.g., scat- 
tering, muon transfer, muonic molecular formation, fu- 
sion), taking into account the dimensions and geometry 
of the apparatus. In particular, it uses the rate of Fig. [y, 
scaled by Cd=0.5, to describe resonant formation in HD. 
The majority of cross sections and rates used in the code 
were taken from [ p| and |12|. However, experimentally 
measured rates for muon transfer (fip — > fit) and nonres- 
onant p/ip molecular formation Q were used (see also 
||). Values for the cross section (J^t+d equal to 0.9 times 
those of Jl^] were found to better represent the fusion 
time distribution observed in the moderation layer [fljl , 
and thus were used in the simulation. These define our 
nominal MC input; the code has been validated through 
comparison with independently written simulation code 
for simplified cases for which a direct comparison can be 
made 

To compare directly the time-of-flight fusion results to 
the simulation, the simulations were scaled by the silicon 
detector solid angle, Qsi = (2.32 ± 0.10) x 10 -2 and nor- 
malized to the fraction of incoming muons which stopped 
in the production layer, S F = (32 ± 2)% @. 

The fusion time-of-flight data for two HD targets are 
shown in Fig. |^. Both the simulation and data are plot- 
ted as fusion events per incident muon N^. Figure ^ (a) 
shows the results for the 16 /ig-cm~ 2 HD target. The sim- 
ulation shows a clear two-peaked structure not evident in 
the data. One peak between 2 and 3 fis corresponds to 
the molecular formation resonances which occur above 
1.0 eV, while the other, between 4 and 5 /is, corresponds 
to the stronger resonances near 0.3 eV. Data from the 
37 /ig-cm -2 HD reaction layer are shown in Fig. ^ (b). 
Here, the fit atoms will undergo more scattering reac- 
tions before molecular formation occurs and hence are 
less sensitive to the molecular formation resonances. This 
is clearly visible in the simulation results, where the two- 
peaked structure is less evident. The agreement between 
simulation and data for this target thickness is quite rea- 
sonable (see below). Good agreement was obtained also 
for the case of a D2 target layer, as reported previously 
1- 

To test sensitivity of our data to resonant molecu- 
lar formation rates, several modifications to the nominal 
MC input were made. Differences from the nominal in- 
put (a) include the following: determining the effect of 
a decreased a^t+d scattering cross section (b); shifting 
the positions of the resonances to higher (c) and lower 
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(d) energies; using an energy independent constant rate 
to model molecular formation (e); and broadening the 
molecular formation resonances (f). Table | summarizes 
the target thickness, input changes to the MC, the fit 
ratio Rf which is the factor by which the data must be 
scaled to provide the best fit to the simulation, and the 
X 2 per degree of freedom (dof=49) for the scaled fit. To 
provide a check of how well the simulation is reproduc- 
ing the absolute intensity observed, the x 2 /dof of fits are 
listed for Rp being held fixed at 1. 

To investigate sensitivity to the positions of the molec- 
ular formation resonances, they were shifted in the simu- 
lation by 0.18 eV (the average standard deviation of the 
peaks) to higher and lower energies. The peaks appear- 
ing in the simulated fusion time distribution shifted as 
expected due to the inverse relationship between the time 
of flight and speed of the fit atom. The results given in 
Table | show increases in x 2 /dof for both the shape and 
intensity tests when compared to the simulation results 
using the nominal predicted positions of the resonances. 
The results are indeed sensitive to the position of the res- 
onances, but shifting them does not improve agreement 
between simulation and data. 

A simulation using a constant molecular formation rate 
of Xdfit = 200 /is -1 , chosen to approximately match the 
intensity of the data, was also compared to data. The 
X 2 /dof for the fit to the 16 fig- cm~ 2 data is 2.99 com- 
pared to 1.87 using the nominal MC input. 

To determine whether the discrepancy in shape could 
be explained by resonance broadening, the theoretical 
resonances were convoluted with a Gaussian of 50 meV 
rms, while preserving the area under each resonance. 
This MC time distribution for the 37 /xg-cm -2 data did 
not reproduce the data as well as the nominal simula- 
tion. For the 16 /ig-cm -2 data, the two-peaked structure 
remained in the MC, with the distribution of both peaks 
quite broad and a factor of 1.6 higher in yield than the 
data for the time range between 3.0 and 6.0 fis. 

To determine the best value for the number (S\) by 
which the resonant molecular formation rate should be 
scaled to agree with the data, simulations were done for 
S\ between 0.50 and 1.40 (S\ = 1.0 corresponds to agree- 
ment with theory) and the results were then fit to the 
data. The x 2 values of these fits versus S\ allow esti- 
mation of the scaling factors for the 16 and 37 /xg-cm~ 2 
targets of, respectively, 0.798 ± 0.060 and 1.040 ± 0.072, 
which are not in agreement within their quoted statis- 
tical errors. As the significant known systematic errors 
in the experiment (such as solid angle and stopping frac- 
tion) are common to both the 16 and 37 /ig-cm~ 2 tar- 
gets, they cannot be used to reconcile the scaling factors, 
pointing to an unknown systematic error or inadequacy 
of the simulation for which an additional systematic un- 
certainty of unknown origin must be determined. Using 
the method described in ]lq ], where the assessed uncer- 
tainty is increased to reflect possible unknown sources 



of error, an uncorrelated systematic uncertainty of 2.38 
times the statistical uncertainty is added in quadrature 
with the correlated systematic uncertainty. The average 
scaling factor for both targets is then found to be 0.897 

±(0.046) stat ± (o.im) syst . 

The scaling factors for the resonant molecular forma- 
tion rates for the targets studied are given in Table [fi|. 
The experimental results confirm the high HD resonant 
rates expected at these energies. However, the fusion 
time spectra structure is not confirmed, in contrast to 
the results for D2 0. The two-peaked structure pre- 
dicted by the Monte Carlo results for the thin HD target 
is not apparent in the data. Several adaptations to the 
resonant molecular formation theory were investigated, 
but this peaked structure persists in the MC. On the 
other hand, using a constant rate to describe molecular 
formation modeled the data less accurately than the res- 
onance model. In the thicker HD target where resonant 
molecular formation competes with scattering of the fit 
atoms and hence energy loss, agreement between theory 
and data is quite good. The theory developed for the 
resonant process does not take into account any effects 
of the crystal lattice or its motion with respect to the 
fit atoms. It has been suggested that effects in the solid 
hydrogen crystal lattice may affect the fit slowing pro- 
cess and thus the reaction yields fTi| . Incorporating such 
effects into the simulation is a logical next step. 
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FIG. 1. Theoretically predicted molecular formation rate 
calculated for 3K, normalized to liquid hydrogen density, for 
the two hyperfine states of 



FIG. 3. Energy deposited in silicon detectors by 3.5 MeV 
fusion alpha particles from the reaction layer, normalized per 
muon (iV M ), for the 16 ^ig-cm -2 HD data. The dashed lines 
show the energy cut imposed for the time-of-fiight analysis. 
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FIG. 4. Time-of-flight spectra for (a) 16 and (b) 37 
/ig-cm -2 HD targets. Data are plotted as solid dots with 
statistical error bars and simulations are shown as histograms 
with the width of the band indicative of the known systematic 
errors of the experiment. 
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FIG. 2. Solid hydrogen layer structure. Fusion can occur 
in either the D2 moderation layer (c), or HD reaction layer 
(d). 
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HD layer = 16 


Mg-cnT 2 






(a) 


nominal (see text) 


0.787(29) 


1.87 


3.23 


(b) 


c^t+d unsealed 


0.839(31) 


1.93 


2.77 


(c) 


+0.18 eV shift 


0.753(52) 


5.86 


7.84 


(d) 


-0.18 eV shift 


0.952(51) 


3.69 


3.69 


(e) 


A dMt =200 /is" 1 


1.03(5) 


2.99 


2.99 


(f) 


convolution 


0.639(21) 


1.63 


9.54 




HD layer = 37 


/ig-cm -2 






(a) 


nominal (see text) 


1.04(4) 


1.27 


1.29 


(b) 


a^t+d unsealed 


1.14(5) 


1.27 


1.47 


(c) 


+0.18 eV shift 


1.06(6) 


2.32 


2.32 


(d) 


-0.18 eV shift 


1.08(5) 


1.40 


1.46 


(e) 


A dMt =200 ms" 1 


0.998(37) 


0.97 


0.97 


(f) 


convolution 


0.751(29) 


1.11 


2.69 



TABLE I. Listed are the changes in Monte Carlo input, 
the fit ratio Rf along with the x 2 /dof for the fit. Fits were 
also done holding the scaling of data to MC fixed at _Rf = 1.00, 
with the resulting x 2 /dof given. 



<~<30 

2 25 

X 

ti 15 
o 

I 10 

\ 5 

+-> 



> ,- 



i l 



2000 3000 4000 5000 

Energy (keV) 



Target 


A dM t Scaling Factor (S\) 


16 ^ig-cm - ^ 


0.798 ± (0.060) stat ± (0.110) corr 


37 jUg-cm" 2 


1.040 ± (0.072) stat ± (O.U8) CO rr 


weighted average 


0.897 ± (0.046) stat ± (0.166) syst 



TABLE II. The values of A^t scaling factors with all 
sources of error included (stat=statistical, corr=correlated 
systematic and sy s£=total systematic error). A value equal 
to 1.0 is predicted theoretically. 
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